Native defects in hybrid C/BN nanostructures 
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First-principles calculations of substitutional defects and vacancies are performed for zigzag-edged 
hybrid C/BN nanosheets and nanotubes which recently have been proposed to exhibit half- metallic 
properties. The formation energies show that defects form preferentially at the interfaces between 
graphene and BN domains rather than in the middle of these domains, and that substitutional 
defects dominate over vacancies. Chemical control can be used to favor localization of defects at C- 
B interfaces (nitrogen-rich environment) or C-N interfaces (nitrogen-poor environment). Although 
large defect concentrations have been considered here (10® cm~^), half-metallic properties can subsist 
when defects are localized at the C-B interface and for negatively charged defects localized at the C- 
N interface, hence the promising magnetic properties theoretically predicted for these zigzag-edged 
nanointerfaces might not be destroyed by point defects if these are conveniently engineered during 
synthesis. 



I. INTRODUCTION 

Hybrid C and BN nanostructures have been studied 
since the first attempts to dope carbon nanotubes with 
boron or nitrogen;- Nanotubes (NT) with C, B, and 
N were theoreticahy predicted^- and experimentally syn- 
thesized in the mid 90s^i^"— but the recent advent of 
graphene has revitalized the field. Hybrid C and BN 
nanosheets are becoming accesible^ ofi'ering a new route 
to enhance graphene's promise, enabling for fine-tuning 
the electronic and optical properties (for example modu- 
lation of the electronic bandgap, chemical reactivity, etc). 

Following the advances in experimental growth of hy- 
brid C-BN nanostructuresji"— the stability of domain- 
separated C and BN nanosheets and nanotubes has been 
studied by first principles simulations in the last cou- 
ple of years Fi2r2i There is general agreement that seg- 
regation is energetically advantageous, because B-N and 
C-C bonds are more stable than C-N and C-B and the 
former are favored by the formation of C and BN do- 
mains. Furthermore, according to molecular dynamics 
simulations, hybrid C-BN armchair nanotubes can be 
spontaneously formed via the connection of zigzag-edged 
BN and graphene nanoribbons (BNNR and GNR) at 
room temperature,^^ and their stability can be competi- 
tive with that of the corresponding pristine BNNTs and 
CNTs.~ These zigzag-edged heterostructures are partic- 
ularly appealing for spintronic applications, as the po- 
larity of the BNNR gives an interfacial dipole, hence an 
effective electric field acting on the graphene ribbonAS It 
has been predicted that half- met allicity can be induced 
in zigzag-GNRs by application of a sufficiently large in- 
plane electric field perpendicular to the edges of the rib- 
bon^ 

Magnetism in zigzag-edged graphene-based nanostruc- 
tures is related to the presence of electronic states that 
are mainly localized at the edgesJ^ These states seem 
to be ubiquitous in hexagonal structures with zigzag 
borders^*, and are not present in armchair edges. Rough- 
ness and defects have been shown to strongly affect the 
electronic properties of zigzag-GNR, inducing a conti- 



nous decrease of the magnetic moments with increas- 
ing concentration of defects until the system becomes 
nonmagnetic for concentrations of one defect per ~10A, 
which might be a typical concentration for real samplesJ^ 
Nevertheless half-metallicity induced by external electric 
fields seems to survive at similar defect concentrations 
with the same critical field strength. At zigzag inter- 
faces between C and BN domains, however, defects can 
change the electronic screening of the edge polarity hence 
affecting the effective electric potential on graphene's rib- 
bon and destroying its intrinsic half-metallicity. It is then 
important to address the role played by point defects at 
the boundary between graphene and BN domains. In 
this work, first principles simulations are performed to 
determine the formation energies and electronic prop- 
erties of intrinsic point defects (vacancies and substitu- 
tional defects) at zigzag interfaces between C and BN 
nanodomains that could exhibit magnetism. It is shown 
that these defects are energetically favored when local- 
ized close to domain interfaces, and their electronic prop- 
erties are determined by the presence of a defect-induced 
electronic state localized at the C-N interface which could 
be occupied or empty depending on the donor/acceptor 
nature of the defect. 

The paper is organized as follows. After describing 
the technicalities in section [Hi we present and discuss 
the findings in IIIII Results of the energetics of neutral 
defects at the interfaces are examined and compared to 
similar defects placed at a graphene-like or BN-like mon- 
odomain. Magnetism at the interfaces and other elec- 
tronic properties are then analysed, and the effect of de- 
fect concentration and charge state discussed. Finally, 
section [IV] summarizes the conclusions. 



II. METHODOLOGY 

Ab initio pseudopotential density functional calcula- 
tions are performed for both planar (2D) superlatticesi^, 
and (n,n) armchair nanotubes (5<n<14) with zigzag 
edged C-BN domains along the tube axis following the 
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FIG. 1: (color-online) Bandstructures for hybrid C-BN nan- 
otubes with different armchair chiralities and zigzag edges 
between C- and BN-domains along the tube axis. Different 
relative compositions of C2,(BN)i_a; are picked to ensure half- 
metallicity. The absence/presence of a band gap for up/down 
spin polarizations is revealed in the solid (red) or dashed 
(blue) lines. 



geometries discussed in the literaturejiiiii^'i^i2£ Troullier- 
Martin type pseudopotentiala^i and numerical atomic or- 
bitals with double-C plus polarization are used to describe 
the electronic valence states within the spin-polarized 
generalized-gradient approximation^ as implemented in 
the SIESTA codej2^ The atomic positions are determined 
with a structural relaxation until the forces are smaller 
than 0.02eV/A. An accurate description of the bound- 
ary electronic states requires a smooth sampling of the 
reciprocal space, and typically Monkhorst-Pack grids of 
at least 1x1x100 is used to sample the Brillouin zone. 

It has been claimed that half-metallicity in C-BN nan- 
otubes can be achieved within a certain ratio of carbon 
content (30%-50%) and is not dependent on the diame- 
ter of the tubei^ This allegation can be misleading. Al- 
though half-metallic ground states can be obtained for 
different chiralities (see FigH]), this only happens for C- 
concentrations under a critical value, that does depend 
on the tube diameter. Analogously to the case of pla- 
nar C-BN heterostructures, the critical parameter is the 
relative width of the C- and BN-strips. -^ To address the 
stability of magnetism to the presence of defects, systems 
that hold half-metallic properties in their pristine (defect- 
free) structure, are considered. In particular, through- 
out this work we take: (i) nanotube with chirality (5,5) 
and 50% C concentration (corresponding to 5 C and 5 
BN zigzag chains), and (ii) planar (5,7) C-BN superlat- 
tice (composed by 5 C and 7 BN zigzag chains). Native 
point defects are introduced near the edges of the C and 
BN domains in a supercell geometry, with at least 30 A 



empty space in-between periodic images. Different con- 
centrations can be modeled by changing the size of the 
supercell along the periodic axis. If not otherwise stated, 
supercells made out of 8 repetitions of the unit cell along 
the tube axis are considered (a total of 160 atoms). 

The formation energy of a defect in charge state q is 
defined in terms of the chemical potentials of the species 
involved in the defect, and the Fermi level, fie, measured 
relative to the top of the valence band Ey : 

nifii + q{fJ,e + Ey) (1) 

i 

where Etot is the energy of the defective supercell con- 
taining Hi atoms of species i (C, B, and N) with chemi- 
cal potential /Xi, at T = (entropic contributions are ne- 
glected). These chemical potentials are specified by a ref- 
erence system that acts as the reservoir of atoms. Typi- 
cally reference values used in the literature are taken from 
graphene, molecular nitrogen, and bulk boron. These, 
however, will only give formation energies under limits 
of high concentrations of C, N, or B. Formally, one can 
use the formation energies of these reference systems as 
thermodyncamic limits to the chemical potentials, which 
must satisfy the following conditions: 

1. stability of the C-BN complex: 

fic + f^B+f^N ^AEf (CBN) (2) 

2. The values that cause precipitation into its con- 
stituents, 

fJ-c < 0, fiB <0, UN <0 

3. The values that cause formation of each domain 

AiB +MW < A£;/(BN) 
Mc < A^^/(C) 

where AEf denotes the generalized formation free en- 
ergy of the corresponding system relative to pure C or 
BN nanotube. Relation 1^ can be used to write one 
of the atomic chemical potentials in terms of the other 
two. Considering that C-doping in BNNT had proved 
elusive^i, whereas BNNT can be synthesized by substitu- 
tional reaction from CNT^^ , in the following the chemical 
potential of C will be set to the value that corresponds to 
clean CNT, meaning that the defect formation energies 
will be given in the C-rich limit, and only fiN will be used 
as free parameter. 

III. RESULTS AND DISCUSSION 

Substitutional carbons on the boron and nitrogen sub- 
lattices (Cb and Cn), boron or nitrogen atoms in carbon 
sites (Be and Nq) and vacancies of each species (Vq, Vb 
and Vn) were considered, as shown in figure O Ghost 
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FIG. 2: (a) Atomic positions considered for the difTerent de- 
fects. Circles, filled squares, and empty squares denote C, 
B and N atoms. Periodic boundary conditions require that 
both C-N and C-B edges are included in the simulation box. 
Formation energies (in eV) for the studied neutral defects as 
a function of Nitrogen chemical potential (/ijv) are shown in 
the lower panel for: (b) the pure CNT, (c) hybrid C-BN nan- 
otubes (solid lines) with 8 repetitions of the unit cell, and 2D 
nanosheet (dashed lines) with 7 repetitions of the unit cell, 
and (d) BN nanotubes with (5,5) chirality. 



orbitals are used to improve the localized atomic orbital 
description at the vacancy sites. Labels N and B will be 
used to denote the edge at which the defect is placed, 
where N designates the interface with C — N bonds, 
whereas B is the edge with G — B bonds. Accordingly 
defects are classified into set N (Cn, Be, V^, Vg , or 
V^), and set B (Cb, Nc, Vg, Vf , or Vg). 

The formation energies for the interfacial (neutral) de- 
fects in (5,5) C-BN nanotube are shown in Figl2]and com- 
pared to the energetics of the same defects in a pristine 
pure-CNT and pure-BNNT of the same chirality. The 
energies of pure C and BN nanotubes are in agreement 
with values published in the literature^^^— , considering 
that the formation energy depends on the chirality, and 
increases with tube's diameter— i^. The cost of remov- 
ing atoms (vacancy formation) is generally larger than 
the energy of atom substitution (disorder), which is in- 
dicative of the high stability of the honeycomb structure. 
First principles calculations have shown that vacancy for- 
mation is energetically more favorable near the edge of 
zigzag BNNRs than in the center of the ribbon22.. A sim- 
ilar effect is observed now, with defects at the edge being 



at least 1 eV less energetic. Interestingly, recent exper- 
iments have shown that carbon substitution (doping) in 
BN nanostructures is favored at the edges.— 

Note that Cn and Be (defects of set N) are energeti- 
cally favorable in a nitrogen-poor environment (B-rich), 
and the same is true for Nq and Cb {set B) in N-rich 
situation, so that the localization of defects at each edge 
of the heterojunction could be chemically controled. Fur- 
thermore, when neutral Cn and Be defects are formed, 
one hole is added to the system, so that these defects 
can act as electron traps. On the other hand, Nc and 
Cb add one extra electron to the system and could pos- 
sibly become positively charged under appropriate con- 
ditions (low electron chemical potential /ie). In pure BN 
nanosheets and nanotubes, Vb is also an acceptor that 
can trap charges to form Vg, whereas Vn has two one- 
electron states in the gap, one occupied and the other 
empty, enabling the formation of both or Vj^, de- 
pending on the value of /ie- 

Figure [3] shows the electronic Density of States (DOS) 
for set N (top) and set B defects. Notice that, even for 
the high defect concentration considered in these simula- 
tions (5-lO^cm^^), the later remain half-metallic, with a 
sharp peak in the DOS at the Fermi level. This peak cor- 
responds to a defect-induced electronic state that is local- 
ized at the N-edge rather than at the defect itself (placed 
at the B-edge). Hence, electrons coming from these donor 
defects are transferred through the C-domain to the re- 
gion with lower electronic potential, as described in ref- 
erence O This can be seen from the DOS projected on 
the carbon orbitals localized at the C-N interface (shaded 
regions in the figure). This localized band has some dis- 
persion along the edge axis (there is slightly more contri- 
bution to the pDOS from C atoms directly opposite to 
the defect site), and a larger unit cell would be needed 
to address its localization length. The exchange splitting 
for this band is ^0.3 eV for all the set B defects, and the 
calculations give a total magnetic moment of '^1 /jb for 
Cb, Nb and Vg, 2 hb for Vg, and ~3 /is for V|. 

On the other hand, the missing electron (extra hole) 
for set N defects has a dramatic effect the correspond- 
ing DOS, that does not show half-metallic properties. 
There is a similar defect-induced electronic state local- 
ized at the C-N interface for Cn, Be and Y^, but now 
it is unoccupied ('^O.l eV above the Fermi level) and 
not spin-polarized (Cn and Be have a total magnetic 
moment of ^0.5 /is). and Vg are slightly different, 
with an empty C-N interfacial state at ~0.4 eV above the 
Fermi level, and a rather localized electronic state on the 
atoms surrounding the vacancy site (shaded light regions 
in the figure) which is spin-polarized (total magnetiza- 
tion of '^1.8 /iB and l/is, respectively. The lowering 
of the Fermi level means a partial depopulation of the 
ttb statCfi^ and a weakening of the magnetic instability 
induced by electron-electron interactionsi^. The reduc- 
tion of the electronic occupation of the ttb state lowers 
the electronic potential at the B-edge, and significant de- 
crease of the induced electric field in the graphene ribbon. 
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FIG. 3: Density of States (DOS) around the Fermi level (at eV) for defects localized at the C-N (top) and C-B (bottom) 
interfaces. Dark shaded regions correspond to the projected DOS (pDOS) on C orbitals localized at the N-interface. Gray 
regions for Vc and Vg are the pDOS on C or N orbitals close to the C or B vacancy, respectively. 



SO that half-metallicity cannot survive>i^ 

Let's focus in the following on those defects that have 
the lower formation energies, i.e. substitutionals. The ef- 
fect of increasing the defect concentration even further is 
shown in figure IH where the DOS of different supercells 
of planar bidimensional superlattices are shown. Small 
differences arise between this planar geometry and the 
previous tubular case, mainly because of the different 
C/BN ratio, but the general results are equivalent. For 
p-type set N defects (Be and Cn) the reduced number 
of electrons shifts the Fermi level towards lower values as 
we move towards higher defect concentrations (smaller 
supercells, from right to left), depopulating the initially 
fully occupied majority spin band, localized at the C-B 
interface (light shaded region), until the full spin polar- 
ization at the Fermi level is destroyed. On the other 
hand, for n-type set B defects (Nc and Cb), the extra 
electrons added to the system shifts the Fermi level to 
higher values as we move towards higher concentrations, 
increasing the occupation of the C-N electronic states, 
and eventually destroying half-metallicity. Notice again 
that defects localized at the C-B interface add extra elec- 
trons that are transferred to the C-N interface, while de- 
fects localized at the C-N interface, which are electron 
deficient tend to remove charge from the C-B interface. 

For the limit of low concentration of defects, the 
semimetallic graphene domain could be considered as 
the electron reservoir so that /Xe is fixed by the posi- 
tion of GNR's Fermi level, and hence only neutral de- 



fect calculations were considered up to now. However 
the donor/acceptor character of the defects has become 
relevant when discussing their electronic properties, re- 
vealing the need to address charged defects. When ad- 
dressing charge states for defects in semiconductors and 
insulatores, the value of /Ze in equation ([T]) ranges be- 
tween zero (the top of the valence band) and the energy 
band gap, and a careful alignment of the reference energy 
has to be considered to compare energetics of neutral and 
charged defects. Here, however, the system is semimetal- 
lic regardless of the charge state, so that the position of 
the Fermi level can be used to align the energetics, and 
fie = can be used to compare formation energies. 

As shown in figure [5l when an extra electron is added 
to the system {q = —1), it localizes at the N-edge state 
(dark peaks) close to the Fermi level, recovering half- 
metallicity for the set N defects (the states become fully 
occupied for set B defects which remain half-metallic). 
The formation energies for these charged defects are re- 
duced by '-^2.3 eV (at /ie = eV). There is a decrease 
in the system's total magnetization for set N (■^O.l (Ib) 
and an increase to 2 fj,B for set B. On the other hand, 
removal of an electron (giving defects with charge state 
q — +l) has negative effects on the electronic properties 
of substitutionals: the lowering of the Fermi level depop- 
ulates the N-edge state for set B defects (although the 
total DOS remains half- metallic), and destroys magne- 
tization for all defects {set remains metallic). In this 
case, the formation energies are increased by ^3.6 eV, 
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FIG. 4: Density of States for substitutional defects as a func- 
tion of defect concentration for planar C/BN superlattices 
(similar results are obtained for nanotube geometries). From 
left to right, supercells with 3, 4 and 7 units of repetition 
along the interfacial axis, corresponding to concentrations of 
13, 10 and 6-10^ cm~^. Dark/light shaded regions correspond 
to projected DOS over C orbitals at the N/B interface. 



clearly destabilizing these positively charged defects. 

Following the above results, a N-rich environment to- 
gether with a high (providing electrons to the system) 
would increase the possibilities for subsistence of half- 
metallicity in these hybrid C-BN structures. B-edge lo- 
calized defects are favored by an N-rich atmosphere and 
even at large defect concentrations (10^ cm~^) conduct- 
ing electrons remain fully spin-polarized. On the other 
hand, under these chemical conditions, there would be 
a lower concentration of defects localized at the N-edge 
(higher formation energies), which although hinder half- 
metallicity when neutral, they do not in their negative 
charge state. It is interesting that similar conditions 
have been suggested to be best for manufacturing hy- 
brid compounds through electron-beam-induced doping 
of BN nanostructures with C— 



IV. CONCLUSIONS 

First principles DFT calculations have shown that in- 
trinsic defects in hybrid C-BN nanostructures are more 
stable at the boundaries between C and BN domains. 



Both bidimensional superlattices of alternating graphene 
and BN nanoribbons, and armchair nanotubes with ax- 
ial domains of C and BN that have zigzag edges have 
been analysed, giving similar results. In these systems, 
the electronic bandstructure close to the Fermi level is 
determined by characteristic zigzag-edge states localized 
at the boundary of the hexagonal lattice, as already dis- 
cussed in the literature. 

The stability of the honeycomb structure results in a 
higher formation energy for vacancies (Vc, Vb, or Vn) 
than for substitutionals. While a nitrogen-poor environ- 
ment favors the formation of defects at C-N interfaces 
(mainly Cn and Be), nitrogen-rich atmospheres favors 
defects localized at the C-B boundary (Cb and Nc), 
which could be classified as electron acceptors and donors 
respectively. Extra electrons provided by the latter are 
transferred to a defect-induced partially occupied elec- 
tronic state localized at the opposite edge (C-N) and the 
system would remain half-metallic, even for the large de- 
fect concentrations considered here. If extra electrons are 
added to the system, they are transferred to this elec- 
tronic state and the defect would be further stabilized 
(in their negatively charged state), whereas electron- 
deficiency drains the electronic level and increases the 
defect formation energy. On the other hand, substitu- 
tional defects localized at the C-N interface, add one hole 
to the system hence lowering the Fermi level and weak- 
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FIG. 5: Density of States for charged (q — —1,0, -fl) substi- 
tutionals. From top to bottom Be, Cn {set N), Cb, and Nc 
{set B). Dark (light) shaded regions correspond to the pDOS 
on C orbitals localized at the N-interface (B-interface) . 
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ening the magnetic instability (through depopulation of 
the T^c-B orbital) resulting in the cancelation of half- 
metallicity. The defect-induced electronic state localized 
at the C-N interface is still there, but unoccupied. Ex- 
tra electrons added to the system would populate this 
level, recovering half-metallicity and lowering the defect 
formation by a couple of eV. 

Finally, the analysis of the magnetic properties reveals 
that half-metallicity (i) could survive defect concentra- 
tions at the interfacial line of up to lO"^ cm"'^, and (ii) 
could be tuned through appropriate electronic doping 
that avoids an excesive charge depletion from the C-B 



interface for high concentrations of p-type substitution- 
als or excesive charge increase at the C-N interface for n- 
type substitutionals. Experimental investigation of mag- 
netism in these hybrid nanostructures is encouraged. 
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